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ABSTRACT 

We present a model for spectrophotometric calibration errors in observations of quasars from the 
third generation of the Sloan Digital Sky Survey (SDSS-III) Baryon Oscillation Spectroscopic Sur¬ 
vey (BOSS) and describe the correction procedure we have developed and applied to this sample. 
Calibration errors are primarily due to atmospheric differential refraction and guiding offsets during 
each exposure. The corrections potentially reduce the systematics for any studies of BOSS quasars, 
including the measurement of baryon acoustic oscillations using the Lyman-a forest. Our model sug¬ 
gests that, on average, the observed quasar flux in BOSS is overestimated by 19% at 3600 A and 
underestimated by ^ 24% at 10,000 A. Our corrections for the entire BOSS quasar sample are publicly 
available. 

Keywords: cosmology: observations, intergalactic medium, quasars: absorption lines, quasars: general 


1. INTRODUCTION 

One of the main goals of the third generation of the 
Sloan Digital Sky Survey (SDSS-III; Eisenstein et al 
2011|) Baryon Oscill ation Spectroscopic Survey (BOSS 


awson et al. 


2013) is to measure the baryon acoustic 


oscillation (BAO) scale in the Lyman-a for^ 
servations of high-redshift {z > 2) quasars ( 

est from ob- 

Busca et al. 

2013 

ISlosar et al.||20131 Kirkby et al.||2013| 

Lbnt-Ribera 

et al.||2014| Delubac et al. 2015|). The transmitted flux 


hydrogen den sity (Gunn & Peterson 1965 Rauch 1998 
Meiksin 2009) along the line of sight that can be nsec. 


to infer th e clustering of the underlying dark matter dis¬ 
tribu tion dCen et al. 1994| Bi et al.||1995| |Zhang et al 


199^ [Hernqnist et al 


nuMim „ 

Bi fc bavidsen|1997 Hui et al.|1997[ Thenns et al.|i998 ^ 

Measurements of the BAO scale constrain the expansion 

history of the universe and can be used to infer t he char¬ 
acteristics of dark energy (Weinberg et al.||2013). 

During the period 2009-14, BUSS observed 294,512 
quasars from a total sky area of 10,400 square degrees. 
As a result, the survey contains the largest sample of 
spectroscopic quasar observations to date and enables 
an unprecedented view into the mul tiple areas of quasa r 
science, e.g. clustering of quasars ([White et al.|2M^, 
quasar luminosity function dPalanqne-Dela bronille et al 
2013 Ross et al.|2013 McGreer et al.|2013| ), and variabil- 
ity properties of broad absorption l ines ~ni quasar spectra 
caused by high-velocity outflows (|Filiz Ak et al.[ 2012 


|1995| Zhang e 
iralda-Escnd^ et al.||19961 


2013[ Hall et al. 2013[ Vikas et al.||2013[ ~ l^Jliz Ak et al 
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Une potential source of systematics for any study of 
BOSS quasars is spectrophotometric calibration errors. 
Calibration error s in BOSS are larger than for SDSS-I 
( York et al.|2QQQ ) due to a design tradeoff that improves 
throughput in the Lyman-a forest of quasar spectra. In 
this work, we describe the dominant source of these er¬ 
rors in BOSS spectra and our procedure for reducing 
them. The miscalibration of BOSS spectra, on average, 
accounts for a ^ 19% excess at 3600 A and a ^ 24% 
decrement at 10,000 A with a smooth transition between 
(see Figure]^. 

This paper is organized as follows. In Sectionwe de¬ 
scribe the SDSS data samples used in this analysis, along 
with the relevant charactistics of the survey design. We 
describe our methods for modeling and correcting the 
spectrophotometric calibration errors present in BOSS 
quasar spectra in Section In Section we evaluate 
our corrections using observations from a BOSS ancillary 
program that contains offset spectrophotometric stan¬ 
dard stars, repeat observations of quasars observed in 
both SDSS-I and BOSS, and a sample of stellar contam¬ 
inants in the BOSS quasar sample. Finally, we conclude 
with a discussion of our findings in Section 

2. DATA SAMPLES 

We us e data from the v5_7_0 BOSS spectroscopic 
pipeline ( Bolton et al.|2Ql^ pr ocessing of the SDSS Data 
Release 12 (Alam et al. 2015|l[^ There are two primary 
classes of targets used in this work: spectrophotometric 
standard stars (main s equence F stars used f or ca libra¬ 
tion) and q uasars; see |D awson et al.| (|2013|) and [Ross 
et al.| (|2012), respectively, for descriptions of the target 
selection for these samples. 

The BOSS double spectrograph mounted at the 
Cassegrain focus of the 2.5-m SDSS telescope (Gunn 
et al. 12006 ) located at Apache Point Observator 3 l^(APO) 
simultaneously records 1000 spectra over a 3^field of 
view. The light from each target is captured by an op- 


^ http://www.sdss.org/drl2/ 
® http://WWW.apo.nmsu.edu/ 
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tical fiber plugged into an aluminum plate at the fo¬ 
cal plane and trans ported to either one of two spectro- 


spc _ 

graphs for analysis ( [Smee et aL|2Q13| ). The circular focal 
plane hole for each 2 arcsecond (120/im) diameter fiber 
is predrilled according to each target’s sky position at a 
nominal observing time. Each BOSS plate includes 20 
spectrophotometric standard star targets which are used 
for spectrophotometric calibration. The calibrations are 
derived by fitting stellar spectrum models to the recorded 
spectra for these targets. Spectrophotometric calibration 
errors in BOSS quasar spectra are primarily due to off¬ 
sets in fiber hole positioning between quasar targets and 
spectrophotometric standard stars that are intentionally 
introduced to improve the signal to noise ratio of the 
Lyman-Of forest region of high-redshift quasars. Addi¬ 
tionally, there are 16 “guide” stars for each plate which 
are each observed with coherent fiber bundles and used 
for guiding the telescope during exposures; see |Dawson| 


et al. (2013) for more details. 

The BOSS DR12 sample contains a total of 487,276 
targets with focal plane offsets, distributed between 2,377 
observations of 2,340 plat es. Of those tygets. the BOSS 
data processing pipeline (Bolton et al. 2012) identified 
284,085 as quasars and 159,886 as stars! These stars 
are often referred to as “failed quasars”, since they were 
targeted as quasa rs due to their photometric similarities 
( Ross et al.[|2012 ). Table summarizes the different tar- 
get samples. 

In addition to the primary BOSS spectrophotometric 
standard stars, a small subset of observations contain 
a second sample of spectrophotometric standard stars 
which were selected and designed to have the same off¬ 
sets in the focal plane as the BOSS quasar targets. This 
sample was collected as part of an anc illary program 
(Dawson et al. 2013 Alam et al. 2015). These stars 
were selected and visually inspected to ensure a uniform 
distribution across the focal plane, in a manner similar 
to normal spectrophotometric standard stars targeted in 
BOSS. The algorithm for photometric selection applied 
to this sample was identical to the algorithm for primary 
spectrophotometric standard stars in BOSS. 

Spectrophotometric standard stars offset in the focal 
plane are identified by bit 20 of the ancillary_target2 
keyword in DR12 data modej^ We will refer to these ob¬ 
jects as “offset standards”. There are 1,770 offset stan¬ 
dard star targets confirmed as stars distributed between 
161 observations. Of those observations, 79 contain at 
least 10 offset standards in total, and 20 contain at least 
10 offset standards per spectrograph. 

As part of this ancillary program, we modified the 
BOSS data processing pipeline to use the offset stan¬ 
dards for spectrophotometric calibration, instead of the 
normal spectrophotometric standard stars. This sample 
provides a vital cross-check of the throughput correction 
model described below. In particular, we use the sam¬ 
ple of 20 plates with at least 10 offset standards in each 
spectrograph for validation tests. In Tablethe column 
labeled Validation indicates the number of targets in the 
validation sample for each of the relevant target samples. 
We refer to the sample of self-calibrated offset standard 
stars as “offset standards”. 

The design wavelength Xi for a target is encoded in 


the LAMBDA_EFF quantity of a plateDesign file and is 
propagated to a target’s corresponding entry in spPlate 
and spAll files in the DR12 data modepq The design 
hour angle ho for each plate is encoded in tEe HA quantity 
of its plateHoles file, but is not propagated to spPlate 
or spAll. We calculate the design altitude uq using the 
design hour angle and the central sky position of a plate. 

The distribution of point-spread function (PSF) full- 
width half-maxima (FWHM) for all observations in 
DR12 is shown in Figure For each observation, we 
summarize the PSF FWHM using an unweighted mean 
calculated from the individual exposures of the obser¬ 
vation. The PSF FWHM is available via the SEEING50 
keyword in spCFrame files for each exposure and is esti¬ 
mated from guide star images. Similarly, we calculate the 
mean hour angle for an observation hobs and the mean 
altitude Uobs of fhe individual exposures. We calculate 
the observing hour angle from the midpoint of the times 
stored in taibeg and taiend keywords in spCFrame files 
for each exposure. We display the distribution of hobs for 
plates relative to their design observing hour angle ho in 
Figure 



Figure 1. Distribution of focal plane PSF FWHM, as indicated 
by the seeingSO keyword, for DR12 observations. A total of 79 
observations are missing this information and omitted here. The 
mean (median) FWHM is 1.54 arcseconds (1.49 arcseconds). 

The BOSS data processing pipeline combines individ¬ 
ual 15-minute exposures, typically about four, into a sin¬ 
gle co-added spectrum. Ideally, the individual exposures 
should be corrected prior to co-addition, however, few 
BOSS analyses currently make use of individual expo¬ 
sures. Instead, we calculate a correction for co-added 
data using the mean PSF FWHM and observing hour 
angle of the exposures for an observation. There were 79 
plates observed during the beginning of the survey which 
do not have PSF FWHM data recorded; for those obser¬ 
vations, we use the sample median as an estimate for the 
PSF FWHM. 

3. METHODS 
3.1. Fiber illumination 

The nominal transformation from sky coordinates 
(a, 6) to focal-plane coordinates r = (x, y) depends on: 


http://www.sdss.org/drl2/algorithms/bitmasks/ 
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Target Sample 

LAMBDA EFF (A) 

OBJTYPE 

CLASS 

A^DR12 

-^validation 

Offset targets 

4000 

- 

- 

487,276 

4,104 

Quasars 

4000 

QSO 

QSO 

284,085 

1,737 

Failed quasars 

4000 

QSO 

STAR 

159,886 

1,049 

Spec, standards 

5400 

SPECTROPHOTO_STD 

STAR 

49,635 

400 

Offset standards 

4000 

- 

STAR 

1,770 

486 


Table 1 

Target samples used in this work. For the DR12 samples, we include any potentially useful spectra, only requring that the fiber was 
plugged (i.e., that the unplugged bit of the zwarning bitmask is not set). For the validation samples, we additionally require that the 
spectra have no known problems (zwarning = 0) and that the spectra are from plates in the validation sample (plate keyword matches 
a plate in the validation set). The offset standards are identified via the QSO_STD bit of the angillary_target2 bitmask. 



Figure 2. Distribution of hobs “ ho for DR12 observations. The 
mean and RMS are 2.5° and 9.4°, respectively. 


• The sky coordinates ((ao,(^o) of fho focal plane ori¬ 
gin r = 0. 

• The wavelength A of incident light. 

• The time of the observation, expressed as the local 
hour angle h for the right ascension ao. 

The last two are due to wavelength-dependent refraction 
through the atmosphere with a time-dependent angle of 
incidence equal to the telescope altitude. We assume 
constant nominal conditions for atmospheric tempera¬ 
ture and pressure for all refraction calculations. The at¬ 
mospheric differential refraction (ADR) angle between 
4000 A and 5400 A is about 0.5 arcseconds (the signifi¬ 
cance of these wavelengths is explained below), which is 
comparable to the 2 arcsecond diameter of BOSS fibers 
and the typical point-spread function (PSF) full-width 
half-maximum (FWHM) during DR12 observations of 
1.5 arcseconds. In Figure we show differential refrac¬ 
tion angles relative to both 4000 A (blue) and 5400 A 
(red) light at various observing altitudes. The magni¬ 
tude of ADR decreases at higher observing altitudes. 

Each BOSS plate is designed for a specific pointing 
(q^o, <^o) at a nominal design hour angle Hq. The nominal 
orientation of the focal-plane coordinate system at hour 
angle h is fixed such that the direction is aligned with 
increasing 6 and -\-x is pointing eastwards. Photons at 
wavelength A observed at time h from a target i at sky 
coordinates (a^, 5i) then have a PSF centered at the focal 
plane position 

r^A, h) = r(A, h; 5i, q^o, (^o) • (1) 



Figure 3. Differential refraction of light in the atmosphere. The 
blue (red) curves show the differential refraction of light relative 
to 4000 A (5400 A). Lines show observing altitudes of 50°, 60°, 
70°, 80°, and 90°. At 90°, the differential refraction is zero in 
both cases. The blue and red vertical dashed lines indicate the 
two relevant design waveslengths, 4000 A and 5400 A, respectively. 
The black vertical solid line indicates the interesection wavelength 
which is essentially independent of observing altitude. 

The fiber hole for each target i is positioned at r(Ai, ho) 
so that the the target’s nominal design wavelength A^ is 
centered on the fiber, and therefore has the maximum 
possible throughput, at time ho- Any departures from 
A = Ai or h = ho will therefore introduce a fiber centering 
offset 

dAA,A,,h) = |rAA,h)-r(A,,ho)| (2) 

and a correspondingly reduced throughput that we cal¬ 
culate below. 

In practice, the majority of fibers are positioned using 
the same central wavelength Af = 5400 A. The exception 
are fibers assigned to high-redshift quasar targets, which 
use a bluer A^ = 4000 A in order to improve the through¬ 
put (and hence also signal-to-noise ratio) in the Lyman-a 
forest region (ALya(l+^) = 4000 A at z = 2.29). Figured 
shows an example of r*(A, h) trajectories (the asteriskea 
quantity is explained in the following section) plotted for 
a single quasar target with A^ = 4000 A, covering a range 
of wavelengths A and hour angles h. The origin of the 
calibration errors that we seek to fix in this paper is this 
difference of design wavelengths. Focal-plane offsets were 
first introduced in the SDSS-III BOSS survey, so quasar 
spectra recorded during SDSS-I do not have these cali¬ 
bration errors. 

The overall system throughput for target i includes a 
geometrical factor A^A, A^, h) that measures the fraction 
of incident light that enters its aperture. This fiber ac- 
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Figure 4. Target centroid position relative to a fixed fiber hole as 
a function of wavelength A and observing hour angle h, as predicted 
for an ideal guiding model. The black circle represents a fiber hole 
(2 arcsecond diameter). The colored lines indicate the target’s 
centroid position for specific wavelengths of incident light (4000 A 
to 7000 A in 500 A intervals) as a function of observing hour angle 
h. The observing hour angle window shown here is ± 3 hours from 
the design hour angle ho of the plate. The spacing between points 
is 1 hour. The fiber hole is positioned such that Aj = 4000 A light 
is centered at the design hour angle ho- 

ceptance fraction depends on the PSF size and shape and 
its centroid offset di{X^h). In the following, we assume 
that all PSFs have circular Gaussian profiles described 
by a standard deviation a that is independent of wave¬ 
length and constant during an exposure, so that 

^i(A, A^, h) = ^Gauss(<^5 <^i(A, h)) (3) 

with 

rDi2 

^Gauss(cr, d) = (rd/cr^) dr , 

Jo 

(4) 

where D = 2 arcseconds is the BOSS fiber diameter and 
Iq is a modified Bessel function of the first kind. We 
show later that our corrections are relatively insensitive 
to these assumptions since they involve fiber-acceptance 
ratios. Figureshows the distribution of PSF sizes mea¬ 
sured for DR12 observations and Figure shows exam¬ 
ples of acceptance fractions calculated under different as¬ 
sumptions. 

3.2. Telescope guiding 

BOSS exposures are typically 15 minutes long, dur¬ 
ing which time an active telescope guiding loop uses 16 
bright stars in the focal plane to maintain the pointing 
towards ((ao,do). Without any additional adjustments, 
the variation of di(A,A^,h) during the exposure would 
noticeably reduce target throughputs relative to their op¬ 
timum values at ho where di{Xi, Xi^ho) = 0. Therefore 
the guide loop also makes small adjustments to the focal 
plane position dr, radial scale s, and rotation 0 during 
each exposure, which modify the transformation from 



Figure 5. Acceptance fraction as a function of offset distance 
for typical PSF sizes. The dotted, solid, and dashed curves re¬ 
spectively correspond to acceptance fractions calculated using PSF 
FWHM values of 1.2, 1.5, and 1.8 arcseconds. The green and black 
curves respectively correspond to Gaussian and Kolmogorov PSF 
shapes. The solid vertical line indicates the BOSS fiber aperture 
radius. 

sky to focal-plane coordinates according to 

r ^ G(s, h) • r-b dr , G{s,0) = {1s) R{0) (5) 

where R is a 2D rotation matrix acting on focal-plane 
coordinates. The corresponding centering offset for fiber 
i is then 

d[{X,Xuh) = 

I G{s{h),0{h)) ■ ri(A, h) + 5r{h) - ho) \ , (6) 

where the prime denotes the effects of the guiding ad¬ 
justments dr(h), s(h), and 0{h). 

It is not feasible to reconstruct the actual history of 
guiding adjustments dr(h), s(h), and 0{h) during an ex¬ 
posure so we instead adopt an ideal guiding model, which 
assumes that the plate center tracks (Go,do) optimally 
and that the guider makes continuous adjustments to 
minimize the chi-square figure of merit 

9 

where the sum is over the 16 guide stars, indexed by g 
and given equal weight, and the resulting adjustments 
optimize throughput at A^ = 5400 A. We denote the tra¬ 
jectories resulting from ideal guiding as r*(A, h) and the 
corresponding centering offsets as d*{X^ Xi^h). Figure]^ 
shows an example of trajectories r^A^, h) (left panel) ana 
r*(A^, h) (right panel) at each target’s design wavelength 
Xi for an entire plate. The middle panel shows trajecto¬ 
ries when only the plate center (gq, ^o) is tracked and no 
additional adjustments are made to track the guide stars. 
We find that, under the ideal guiding assumption, the 
translation, scale, and rotation degrees of freedom com¬ 
pletely remove the monopole and dipole contributions 
to centering offsets for targets with Xi = Xg = 5400 A 
(shown in red), but lead to larger offsets with some resid¬ 
ual monopole and dipole for the smaller set of quasar 
targets with Xi = 4000 A (shown in blue). 
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Figure 6. The ideal guiding model for an example plate (PLATE = 6114). The black circle indicates the plate’s 3° field of view. The 
colored lines represent a target’s Xi light path spanning a range of ± 3 hours from the design hour angle of the plate. Red (blue) lines 
indicate Aj = 5400 A (4000 A) targets. The fiber centroid offset distances have been exagerated relative to the plate dimensions with the 
scale indicated in the bottom left of each panel. 


3.3. Spectrophotometric calibration 

The mean rate of detected photons from source i at 
wavelength A and hour angle h is related to the source’s 
true spectral energy distribution (SED) fi{X) by a fiber- 
dependent calibration factor C 

cJti ■ 


standard-star calibrations Cc to other targets in the fo¬ 
cal plane does not account for differences in each target ^s 
design wavelength A^. Since the calibration stars have 
Ac = 5400 A, this only affects the quasar targets, which 
are designed with Xi = 4000 A. Under the ideal guid¬ 
ing assumption, the SED mis-calibration of targets with 
A^ 7 ^ Ac is the ratio of acceptance fractions at the actual 
{Xi) and assumed (Ac) design wavelengths 


which we split 

Ci{X,Xi,h) = Al{X,Xi,h)Bi{X,h) (9) 

into the fiber-acceptance fraction A*{X,Xi,h) assuming 
ideal guiding (denoted by the asterisk), and a term 
Bi that includes all other contributions to the signal 
throughput but is independent of a target’s design wave¬ 
length Xi. The quantities that we actually measure from 
an exposure are the integrated photon counts 


^ A* (A, Ac/tobs) 

“ A*{X,Xi,hobs) ’ 


(13) 


and our improved calibration consists of scaling the 
pipeline SED estimates by this factor 


fi{X) Rf (A) fi{X) . (14) 

To summarize, the corrections i?*(A) that we implement 
depend on: 


/•^obs+A/i/2 7 

nAA,A,,hobs)= / ^{\Xi,h)dh (10) 

Jhobs-Ah/2 

dn ■ 

-^{X,Xi,hobs)Ah, ( 11 ) 

where hobs is the exposure midpoint hour angle and 
Ah = 3.75° (15 minutes) is the exposure duration. The 
variation of the integrand is sufficiently small over Ah 
that we use the trapezoidal approximation to the inte¬ 
gral in the following. Eigure shows the distribution of 
^obs — ^0 for all DR12 observations. 

The 20 spectrophotometric calibration targets on each 
plate, indexed by c, are chosen so that their true SED 
/c(A) can be directly estimated as /c(A) using a stellar- 
model fit to their observed photon counts nc(A, Ac, hobs)- 
We use these SED estimates to provide estimates of the 
calibration factors Cc(A, Ac, hobs) at 20 locations across 
each plate, and then interpolate these factors to the lo¬ 
cations of every fiber, C^A, A^, hobs)- We then estimate 
the SED of all targets on a plate as 

/j(A) = Ci{X, Xi, hohs)~^ ni{X, Xi, /iobs)/A/i. (12) 

A key point of this paper is that the method used 
by the BOSS data processing pipeline to interpolate the 


• The actual fiber centering offsets d* (A, A^, hobs) and 
the assumed offsets d* (A, Ac, hobs) at the exposure 
midpoint hobs (both under ideal guiding condi¬ 
tions). 

• The design hour angle ho, which determines the 
locations of each fiber on the focal plane. 


The PSE size a. 


Eigure shows scatter plots of d*(A, A^, hobs) versus 
d*(A, AcT^obs) for all offset targets in DR12 observations 
(hobs and ho), for six different wavelengths A, with con¬ 
tours of Rf (A) superimposed for three different PSE sizes 
(7. Additionally, we show contours of R*(A) for a nominal 
PSE siz e assuming a Ko lmogorov atmospheric turbulence 
model ( Tat arskiill 1961 ) which includes a wavelength de- 
pendent sca ling of the nominal PSE EWHM (x A“^/^ 

(Eried 1966[ ) . 

In L igure we show the central 68% and 95% quan¬ 
tiles of R*(^ for all 487,276 offset targets in DRIT In 
general, the corrections R*(A) are smooth, monotonic 
functions that cross unity near 4500 A(as expected, see 
Eigure^. Comparing with the subset of Rf{X) where 
l^obs ~^o| < 1-25° (5 minutes), approximately 10% of 
all observations, we see that the vertical spread at the 
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A = 3600 A 
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A = 4500 A 
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A = 4000 A 
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A = 5400 A 
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Figure 7. 2D histrogram of d*(A, A^,/lobs) and d*(A, Ac,/lobs) focal plane position offsets relative to the fiber hole center for all offset 
targets. The color of each 2D bin in a panel corresponds to the number of entries as indicated by the panel’s adjacent vertical color bar. 
The green and black contours represent levels of i^*(A) calculated assuming Gaussian and Kolmogorov PSF shapes, respectively. In the 
panel labeled A = 5400 A the additional dotted, solid, and dashed contours represent levels of R*{X) for PSF FWHM of 1.2, 1.5, and 1.8 
arcseconds, respectively. 
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crossover A is caused by plates observed away from their 
design hour angle. On average, our model suggests that 
the observed quasar flux is overestimated by ^ 19% at 
and underestimated by ^ 24% near 3600 A and 10,000 A, 
respectively. 


3.4. Correction Parametrization 

In order to summarize the general properties of our 
correction model for all DR12 observations, we fit a 
parametrized model 


1 + Cl log A/Ao + C2 (log A/Ao)" 


(15) 


to each throughput correction prediction i?*(A). Fig¬ 
ures [OpT] show the median best fit parameter per plate 
and its relationship to either Uobs or (hobs — ho). The 
points in each figure are colored by (aobs~<^o) for Ihe cor¬ 
responding observation. The Aq parameter, which corre¬ 
sponds to the crossover point of the multiplicative cor¬ 
rection, is strongly correlated with (aobs as seen in 

Figure The value of ci sets the overall size of the cor¬ 
rection and is strongly correlated with the altitude of the 
observation, as seen in Figure Fignre pT] shows a sim¬ 

ilar relationship betwen the magnitude C 2 and observing 
altitude. Observations at higher altitude tend to have 
smaller ci and C 2 values, which corresponds to a smaller 
correction. For reference, we also display the relationship 
between (hobs — ho) and (aobs — ^^o) in Figure 12 


4. VALIDATION 


In order to validate our model, we use observations 
that have been calibrated with offset standards, as de¬ 
scribed previously. In Figure [T^ we display the median 
flux ratio of DR12 failed quasar spectra, with (blue) cor¬ 
rections and without (red) corrections, compared to their 
corresponding spectra from the “validation” reduction. 
The corrections remove a wavelength-dependent system¬ 
atic bias in the distribution of BOSS quasar spectra. The 
vertical spread in each distribution indicates that there 
are other factors affecting the calibration but that they 
are independent of the effect described in this work. Fig¬ 
ure displays a similar comparison between a set of 
quasar targets with spectroscopic observations in both 
BOSS and SDSS-I. In this case, the SDSS-I observation 
is used for reference since focal plane offsets were not in¬ 
troduced prior to BOSS. We note a significant drop off 


below 5000 Abut do not necessarily expect agreement 
between SDSS-I and BOSS because of changes to the 
spectrographs and the intrinsic variability of quasars on 
these timescales. 

Following the asses sment of spectr o photo metric cali- 

20131), we calcu- 


Dawson et al. 


bration performed in ^ _ _ 

late synthetic AB magnitudes ( |Oke Gnnn]|l983 ) using 
SDSS filter-curve- weighteq^ integration of spectra for 
BOSS targets and compareBiem to t heir targets’ corre 


spending SPSS imag ing magnitudes ( |Gnnn et al. 1998 


Fukugita et al. 1996). We apply AB magnitude correc 


tions for SDSS imaging magnitudes {jnAB—'^SDSs) used 
by the BOSS data processing pipeline for the ^-, r-, and 
i-band of 0.036, 0.015, and 0.013. The BOSS spectropho¬ 
tometric calibration of standard stars yields {g — r) col¬ 
ors that are on average 0.022 ± 0.003 magnitudes redder 


^^ http://WWW.sdss.org/instruments/camera/#Filters 


compared to their SDSS imaging magnitudes. Note that 
we calculate the sample dispersion using the normalized 
median absolute deviation. For the offset standards, the 
{g — r) colors are 0.117 ± 0.004 magnitudes bluer before 
applying the correction and 0.006 ± 0.004 magnitudes 
redder after. Using spectra from the validation pipeline 
yields {g — r) colors that are 0.010 ± 0.003 magnitudes 
redder. Fignre [T5| displays the distribution of color resid¬ 
uals A{g — r) between SDSS imaging and BOSS synthetic 
magnitudes for the offset standards. The red and blue 
distributions correspond to synthetic magnitudes calcu¬ 
lated from the BOSS DR12 spectra before and after ap¬ 
plying the corrections, respectively. The black dashed 
line corresponds to synthetic magnitudes calculated us¬ 
ing spectra from the validation pipeline that uses offset 
standard stars for spectrophotometric calibration. 

A complete listing of magnitude and color residuals for 
the relevant data samples is given in Table Our cor¬ 
rections generally remove a significant bias in the mean 
and reduce the dispersion of the magnitude and color 
residuals. 


5. DISCUSSION 

The corrections described here provide a significant im¬ 
provement over the uncorrected quasar spectra in DR12 
and are publicly availablj^ However, there are limi¬ 
tations to attempting to correct the quasar spectropho¬ 
tometric calibration as a post-processing step after the 
data processing pipeline, since many simplifying assump¬ 
tions are required. In order to achieve better results, it 
is necessary and desirable to use as much of the exist¬ 
ing pipeline calibration algorithms as possible but with 
spectrophotometric standards observed with the same 
focal-plane offsets as the quasar targets. Short of in¬ 
cluding spectrophotometric standards with focal plane 
offsets, the throughput corrections described here can be 
applied as a preprocessing step of the calibration pro¬ 
cess. Other improvements over the method described 
here include using a more realistic description for the 
PSF, using the actual guiding adjustments for determin¬ 
ing centroid offsets, and applying corrections to individ¬ 
ual exposures prior to co-addition. This is especially rel¬ 
evant for the extended Baryon Oscillation Spectroscopic 
Survey (eBOSS)[^of SDSS-IV, which uses similar focal 
plane offsets for quasar targets. We are currently work¬ 
ing with the eBOSS data processing team to implement 
the throughput corrections described for individual ex¬ 
posures into the data processing pipeline. Th e proposed 
Dark Energy Spectroscopic Instrument (DESI; Levi et al. 
2013) includes an atmospheric dispersion compensator m 
its design so the effects of ADR should be much reduced 
for that survey. 


Eunding for SDSS-III has been provided by the Alfred 
P. Sloan Eoundation, the Participating Institutions, the 
National Science Eoundation, and the U.S. Department 
of Energy Office of Science. The SDSS-III web site is 
http://www.sdss3.org/, 

SDSS-III is managed by the Astrophysical Research 
Consortium for the Participating Institutions of the 

http://darkmatter.ps.uci.edu/tpcorr/ 
http:// WWW .sdss.org/surveys/eboss/ 
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Figure 8. Throughput correction summary. The shaded regions correspond to the the central 68% and 95% quantiles of the predicted 
correction at each wavelength. The black solid and dashed lines similarily correspond to the central 68% and 95% quantiles of R* (A) for 
observations where |hobs — ^o| < 1-25° (5 minutes), approximately 10% of all observations. 


Data Sample 

# Targets 

Ma 

Mean 

-r) 

Disp. 

A(r 

Mean 

Disp. 

Mean Disp. 

Ar 

Mean 

Disp. 

Ai 

Mean 

Disp. 

Spec, standards 

400 

0.022 

0.057 

0.014 

0.032 

0.032 

0.057 

0.011 

0.052 

-0.004 

0.060 

Offset standards 

486 

-0.117 

0.098 

-0.053 

0.068 

0.084 

0.122 

0.201 

0.159 

0.255 

0.214 

Corr. offset standards 

486 

0.006 

0.085 

0.004 

0.062 

0.047 

0.124 

0.041 

0.163 

0.037 

0.198 

Spec, offset standards 

486 

0.010 

0.077 

0.008 

0.058 

0.054 

0.086 

0.044 

0.097 

0.036 

0.127 

Failed quasars 

1,049 

-0.067 

0.132 

-0.065 

0.127 

0.071 

0.158 

0.138 

0.187 

0.203 

0.245 

Corr failed quasars 

1,049 

0.058 

0.135 

-0.006 

0.116 

0.029 

0.157 

-0.028 

0.186 

-0.023 

0.238 


Table 2 

Band-pass filter magnitude and color differences between SDSS imaging and synthetic AB magnitudes calculated from BOSS spectra. The 
sample mean and dispersion are listed for various data samples. Dispersions are calculated as the normalized median absolute deviation. 



Figure 9. Relationship between median fit parameter Aq and 
hour angle difference (hobs — ho) for all 2377 DR12 plates with 
offset targets. Points are colored by the corresponding altitude 
difference (o-obs ~ <^ 0 ) for the observation. 


Figure 10. Relationship between median fit parameter ci and 
observing altitude o-obs for all 2377 DR12 plates with offset targets. 
Points are colored by the corresponding altitude difference (aobs “ 
ao) for the observation. 
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Observed Wavelength A (A) 

Figure 13. Median flux ratios of failed quasar spectra from the 20 validation plates. The blue (red) curve shows the median flux ratio 
between corrected (uncorrected) BOSS spectra and their corresponding “validation” calibrated spectra. The shaded region corresponds to 
the 68% level of the distribution in each wavelength bin. 
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Figure 14. Median flux ratios of quasar spectra that were observed in both SDSS-I and BOSS. The blue (red) curve shows the median 
flux ratio between corrected (uncorrected) BOSS spectra and their corresponding SDSS-I spectra. 



A(p-r) 

Figure 15. Residual {g — r) color distribution between SDSS 
imaging and BOSS synthetic magnitudes for a sample of offset 
standard stars. The red and blue distributions correspond to syn¬ 
thetic magnitudes calculated from the BOSS DR12 spectra before 
and after the correction, respectively. The black dashed line corre¬ 
sponds to synthetic magnitudes calculated from the modifled data 
processing pipeline that uses offset standard stars for spectropho- 
tometric calibration. 





















































































































